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Abstract 

Theory of nonlinear resonance, including stochastic one, is developed on the basis 
of the statistical field theory and using variables action-angle. Explicit expressions 
of action, proper frequency and nonlinearity parameter as functions of the system 
energy and the external signal frequency are found for the cases of nonlinear pen- 
dulum and double well potential. 
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1 Introduction 



Since the late 1980s, a variety of theoretical and experimental papers appears devoting to 
the stochastic resonance study and discovering new applications in different fields of science 
and engineering [1]. Nowadays, the stochastic resonance is a well established phenomenon 
displaying in a bistabile system simultaneously driven by noise and a periodic signal. There 
is appeared to be an optimal noise level at which the system will exhibit almost periodic 
transitions from one state to other with the frequency of the coherent signal. The enhance- 
ment of a weak input signal has been suggested to characterize by the signal-to-noise ratio 
(SNR) that takes a maximum value at an optimal noise level, i. e., a behaviour which is 
reminiscent of a usual resonance phenomenon. 

The observed enhancement is not due to the matching of two frequencies, but rather to a 
cooperative effect of the coherent signal and the noise. This peculiarity is reflected formally 
as follows. The kinematic condition 
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IXTk = 



(1) 



relates the frequency of the coherent signal f2 to the Kramers rate rx, which is the function 
of the noise level type of temperature T: 



Here, prime denotes the derivation with respect to a generalized coordinate x, F(x) is a 
related dependence of the effective energy, which has a minimum at the point Xq, a maxi- 
mum at Xb and a barrier of the height AF = F(xb) — F(xq); m is a particle mass, 7 is a 
kinetic coefficient. Physically, the kinematic condition means that for every frequency Q, the 
stochastic system could picks out suitable noise level T. However, the value SNR [1] 



being found on the basis of the statistical theory, displays the stochastic resonance maximum 
at the temperature 



which value is fixed by the energy barrier AF, but not the external signal frequency fl. 

It is easily to see that the reason of this contradiction is using the linear approach for the 
SNR determination. In this paper, we develop the theory of nonlinear resonance, including 
the stochastic one. The main ingredients of our approach are: (i) the noise accounting by 
means of introduction of a generalized momentum and (ii) passage to the variables action- 
angle, being usual at studying nonlinear phenomena [2]. Section 2 is devoted to the first of 
these approaches on the basis of the statistical field theory [3]. Section 3 contains details of 
examination of nonlinear resonance in terms of variables action-angle. The latters are tested 
on the example of the simplest model of harmonic oscillator. Then, explicit expressions of 
action, proper frequency and nonlinearity parameter as functions of the system energy are 
found for the cases of nonlinear pendulum, double well potential and nonlinear pendulum 
under constant external field. Assuming the proper frequency be identical to the frequency 
of external signal, we examine resonance conditions of the nonlinear pendulum in Section 4 
and the case of the stochastic resonance in Section 5. Conclusion in Section 6 shows that the 
above mentioned contradiction is resolved because the nonlinear resonance condition fixes 
the proper frequency as function of the external one, but not the noise level, as in the linear 
case. 




(2) 



SNRocT- 2 r K (T), 



(3) 



T m = l -AF, 



(4) 



2 



2 Basic equations 



Let us study a hydrodynamic mode amplitude x(r,t) which space-time dependence is deter- 
mined by the Langevin equation 

x(r,t) = 1 f(x,t) + ((r,t), (5) 

which can be interpreted in the Ito sense. Here the dot indicates the differentiation with 
respect to time t, r is coordinate, 7 is a kinetic coefficient, f(x,t) is a force conjugated to 
the stochastic variable x, and ((r,t) is a stochastic term having the form of white noise: 

(C(r, t)) = 0, (C(r, t)C(0, 0)) = T5(r)5(t), (6) 

where the angle brackets denote averaging and T is a noise intensity type of temperature. 2 
For the system under consideration, the total force 

f(x,t)=f (x) + f ext (t) (7) 
consists of usual internal term 



F(x) + ^|Vx| 2 



dr, (8) 



where F(x) is a system potential per unit volume, (3 > is a constant, V = d/dr, and 
external harmonic term 



f ext (t) = Acos(Qt + ip) (9) 

being determined by an amplitude A, a frequency Q, and an initial phase ip. It is now conve- 
nient to go over to dimensionless quantities by referring the coordinate r to a characteristic 
spacing a, the time t and the inverse frequency Q" 1 to the scale a 3 /T'~f, the amplitude A, 
the internal force / , and the quantity F to T/a 3 , and the fluctuation ( to T^/a?. In this 
case, Eq.(5) can be written as follows: 



x = [V 2 x + f(x, t)} + C(t), f{x, t) = f (x) + f ext {t); 



2 The accepted choice of the coefficient in correlator (6) makes it possible to preserve the standard 
form of the Focker-Planck equation in which the diffusion term appears with the coefficient 1/2 
(accordingly, the Gaussian distribution has the standard dispersion). The doubling of the coefficient 
leads to the standard Onsager form of the diffusion component of the flux probability. 
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(10) 

f = -dF/dx, f ext (t) = Acos(Qt + ip). 

The range of applicability of the Ginzburg-Landau approximation (8) is determined by the 
condition according to which the scale a is much smaller than the correlation length £ = 
^1/2 \Q 2 F/Q X 2 I^q 2 [4]. Averaging Eq.(5) and disregarding correlations, we obtain the Landau- 
Khalatnikov equation for the order parameter (x(r,t)). 

The standard field scheme [3] is based on the investigation of the generating functional 
corresponding to the stochastic equation (10). It is a functional Laplace transform 

Z {u(r, t)} = J Z {x(r, t)} exp (^J uxdrdt^j Dx(r, t) (11) 

for the generalized partition function 

Z {x(r, t)} = ( J] 5 {x(r, t) - V 2 x(r, t) - /(r, t) - C(r, t)} det <K ' " ' ] 

\(r,t) 



5x(r, t) 



(12) 



Here, the argument of the (^-function reduces to the Langevin equation (10), and the deter- 
minant, providing the passage from the continual integration over C( r ^) to x(r,t), is equal 
to the unity within the Ito calculus. 

In the framework of the standard approach [3], the n-fold variation of the functional (11) 
with respect to the auxiliary field u(r,t) allows one to find the n-th order correlator for the 
hydrodynamic mode amplitude x(r,t) and to construct the perturbation theory. However, 
we shall proceed from expression (12) for the conjugated functional Z{x(r,t)}, variation of 
which leads to the most probable realization of the stochastic field x(r,t). Obviously, in the 
framework of the mean-field approximation the functional (12) reduces to the dependence 
Z{{x(r,t))}, which corresponds to the Landau free energy F{(x(r,t))} = — T In Z{(x(r, t))} 
[4]- 

Passing to the consideration of the functional (12), we represent the 5-function in the integral 
form 

ioo 

5{x(r,t)}= J exp (- J pxdrdt^J Dp. (13) 



Then, averaging over the noise ( with using the Gauss distribution 

P {C}ocexpj-i|c 2 M)drdt}, (14) 
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which corresponds to the condition (6), and taking into account Eq.(13), we reduce the 
functional (12) to the standard form 

Z{x(r,t)} = J P{x(r,t),p(T,t)}Dp, P = e' s . (15) 

Here, probability distribution P{x,p} is given by action S — J Cdrdt, where Lagrangian is 

C = p (± - V 2 x - /) - p 2 /2. (16) 

Further, we use Euler equations 

dC d dC „ dC „ 2 dC &R r , 

where dissipative function is 72. = ^i 2 . As a result, equations for the most probable realiza- 
tions of the stochastic fields x(r,t), p(r,t) take the form: 

x = (V 2 x + /) + p, (18) 
p = -V 2 p - f'p - x, (19) 

where prime stands for derivation with respect to variable x. Comparison of (18) with the 
stochastic equation (10), having the same form, shows that the field p(r,t) is the most 
probable value of the fluctuations C( r !^) °f the conjugate force. Differentiating Eq.(18) with 
respect to the time and inserting result into Eq.(19), we obtain the equation of motion as 
follows: 



x + (l + f') x = 2f'V 2 x + / + //', (20) 

where only terms of the lowest order of spatiotemporal derivations are kept. 

So, usage of above field theory allows us to pass from the stochastic equation of motion (5) 
of the first order of differentiation to the equivalent system of two equations (18), (19) of the 
same order, or to the single equation (20) having the second differentiation order. Further, 
we need in using Hamiltonian 7i = px — C that depends on the field variable x and the 
conjugate momentum p. According to Eq.(16), Hamiltonian can be written in the form 

Ti(x,p;t) = H (x,p) + Hi(p,t); 

(21) 

Ho = -VxVp + \p 2 + pf , U x = Apcos{Qt + <p). 
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It is easy to see that these expressions, being inserted into dissipative Hamilton equations 



. an „on . fan „&h\ an 

£ = ^-V— -, P = ~ [ t; V——— I — , (22) 



dp <9Vp' \ dx dVx J dx ' 

arrive at the equations of motion (18), (19). 

It is very important to take into account further that the most probable amplitude p of fluc- 
tuations varies near the magnitude —fo, so that we ought to pass to oscillating momentum: 
p+fo — * P- Moreover, as is seen from the Fourier transformation (13), the momentum p takes 
imagine magnitudes, so that the power f has to consider as imagine one, as well, and the 
sign in front of the last term of Tio in Eqs.(21) must be reversed. Q As a result, Hamiltonian 
(21) takes the final form 

H{x,p;t) = H {x,p) +H\{p,t)\ 

(23) 

= \v 2 + \fl ^1 = -Af cos(Qt + <p), 
where gradient terms are suppressed for brevity. 



3 Nonlinear resonance in terms of variables action-angle 



To analyze the system of equations (18), (19), it is convenient to use the phase portrait 
method. However, in our case such a portrait flows in the course of the time due to appearance 
of the time-dependent external force (9). To avoid this time- variation in phase portrait we 
need in passage from above used variables x, p to new ones - action I and angle •& defined 
as follows [2]: 



Ztt J 

(24) 

o 

where the shorted action S(x, I) plays a role of the generalize function. The convenience of 
the so-introduced variables is that the zero Hamiltonian 7i in Eqs. (23) does not depend on 

3 Apparently, such a reverse corresponds to the passage from above used Euclidean field theory to 
the usual one. 
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the angle so that the corresponding phase portrait is stable in the course of the time. The 
equations of motion for the generalized coordinate i? and the conjugate momentum I read 
(cf. Eqs.(22)) 

tf-^ (25) 

* ~ or L ~ (25j 



where non-homogeneity and dissipation effects are suppressed. According to the second of 
these equations the action I is a constant if an external perturbation is absent. 



3.1 Harmonic oscillator 



To recall a physical meaning of the variables i?, I introduced, let us consider firstly the 
simplest case of harmonic oscillator. In this case, the internal power in Eq.(23) 



fo = -uj x (26) 

is linear to be fixed by a proper frequency c^o- Then, the Hamilton equations (22) arrive at 
the equation of damping oscillation under external power: 



x + x + uJqX = — Aujq cos(f2t + (p). (27) 

This equation differs crutially from Eq.(20) because of the former corresponds to the mo- 
mentum origin p — 0, whereas the latter - to p — —fo- According to Eq.(27), dissipation 
shifts resonance frequency from proper magnitude ojq to value 

w = ^ - 2- 2 , (28) 
whereas a maximum real part of characteristic relation x/A relates to frequency 



Umax = \Ju$ ~ 2-K (29) 

Such a character of the dissipation influence keeps at accounting for anharmonicity effects if 
under parameter c^o one means a proper frequency iv(7i) of nonlinear oscillations depending 
on the system energy. 

To demonstrate advantages of using variables action-angle, let us calculate now their mag- 
nitudes at condition that external power in Hamiltonian (23) is switched off. Then, the first 
Eqs.(24), (25) give immediately 
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I = —, •& = U) . 

Respectively, the shorted action and angle take the form: 



(30) 



S = I 



x \ x 

arcsm | — ) H 

x J x 



'1- 



X 



x 



f x \ 2 2/ 
v = arcsm — , x a = — . 

\x J u 



(31) 



The last of Eqs.(30) gives the usual relation between the angle and the time 



■d = u t + -&o, 
which using arrives at the harmonic laws of motion 



(32) 



x = (2H/ujl) 1/2 sm(uj t + 7? ), 
p = (2H) 1/2 cosm(iu t + 7? ). 



(33) 



3.2 Nonlinear pendulum 



The simplest model of system with possibility of the barrier overcoming is known to be the 
nonlinear pendulum (in this Subsection, we consider the pendulum without a friction and 
an external perturbation). Here, Hamiltonian takes the form 



H = ^p 2 + 2u 2 sm 2 (x/2), 



(34) 



corresponding to the power f = — 2u> sin(a;/2) in Eqs.(23). Hamilton equations (22) arrive 
at the system 



x = p, p = —u sin x. 
Combination of these equations gives the nonlinear one: 



(35) 



x + u sinx = 0. 



(36) 
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This equation is non-solvable in analytical form and we ought to use the phase portrait 
method. The form of this portrait, following from Eqs.(35), shows that the system behaviour 
is governed by energy TC with respect to the critical value H c = 2cJq. At condition H < H c , 
the system moves finitely, whereas with overcoming the critical energy Tt c it passes to infinite 
motion. Let us describe such a behaviour quantitatively. 

In this line, the simplest topic is the solution corresponding to separatrix, for which the 
energy is critical one: H = 7i c . In such a case, the definition (34) gives the separatrix form 
as follows: 



p = ±2w cos(a;/2). (37) 
Then, the first of Eqs.(35) arrives at the separatrix law of motion 



x = 4arctanexp(±a;o£) — tt, (38) 

where the different signs correspond to upper and lower branches of the separatrix. This 
dependence can be written in much more elegant form cos(x/2) = [cosh(±u;o£)] -1 , insertion 
of which into Eq. (37) arrives at the famous soliton dependence 



cosh(u;o£) 

where choice of signs corresponds to solitons moving to right or left sides. In accordance 
with the motion laws (38), (39), the system behaviour on the separatrix (37) is as follows: 
at time t = — oo the system is located in a saddles S T , where the coordinate x = =F7r 
and the momentum p — 0. In the course of the time within the domain — oo < t < oo, 
the former arises monotonously from — 7r to ir, whereas the latter increases at t < and 
decreases at t > 0. It is characteristically that coordinate variation and finite magnitudes of 
the momentum take place within the domain At ~ u^ 1 located near the time t — 0. The 
forms of the corresponding kink x(t) and soliton p(t) are depicted in Fig. 2. 

General solution of Eqs.(35) can be obtained with using the variables action-angle defined 
by Eqs.(24). It is convenient to introduce a parameter 



^A-,^^ n c ^l (40) 

taking the magnitude k — 1 at critical energy H = 7i c , and a new variable £ defined by 
equalities 
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k 1 sin(x/2) 
sin(x/2) 



at 
at 



K < 1, 
K > 1. 



(41) 



Then, the generating function is expressed in terms of the incomplete Jacobian elliptic inte- 
grals F(£,k), E(£,k) of the first and second orders [5] as follows: 



S(x, I) = 4cj 



; [^(e, «)-(i- « 2 )F(e,«)] 



V«) 



at k < 1, 
at k > 1. 



(42) 



Differentiation of these equalities with respect to J arrives at expressions for the angle $ 
that generalyzes the last equality (31) (we supress these expressions because of its very 
complicative form). 

Fortunately, formulas for action I = 4S'(£ = tt/2)/2tc follow from Eqs.(42) immediately and 
are expressed by means of the complete Jacobian elliptic integrals K(k) = F(£ = 7r/2, «), 
-E(k) = E(£ = tt/2,k). Taking into account corresponding dependences shown in Fig. 3, we 
shall need further in using asymptotics of these integrals [5] 



K(k) 



I (l + f) 
In 



at k C 1, 
at 1-k 2 «1; 



(43) 



E(k) 



1(1-*) 



1 + 



1-k" 



In 



at /t<l, 
at 1 - k 2 < 1. 



(44) 



Resulting dependence /(7i) depicted in Fig. 4 shows monotonic increase from I — at 
7i = to infinity with the logarithmical inflection at the critical energy 7i c . This behaviour 
is characterized by the following asimptotics: 



2ujq < 



Mi- 1 



H/H c 



-H/Hc 



In 



16 



1-H/HcJ 

tt V 1 ^ 4 in n/Hc-i) 
2, /-^ f 1 — 2^ 



at H« 7Y C , 
at < H c - Tt Tt c , 

at < H - H c H c , 

at H» 7i c . 



(45) 



Accounting the properties of the elliptic integrals [5], we obtain for the proper frequency 
c<j = d determined by the first equality (25) in the following form: 



7T ^0 f 1 

2K(k) U 



at /c < 1, 
at /c > 1, 



(46) 
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where 



"-{«-« at < 47 » 

As is seen in Fig. 4, the proper frequency falls down from the bare magnitude u at the 
minimal energy 7i = to zero at Ji — 7i c and then, after infinitely sharp cusp, the value u 
increases monotonously. According to Eqs.(43), such a behaviour is presented by asimptotics: 

f 1 ~WT C at H«H C , 

lu sa cj < 7T (in i^/^i )' 1 at |W - W c | < W c , (48) 

On the other hand, definitions (34), (41) arrive at time-dependencies of the momentum: 

f cos£ = cn(t, re) at H < H c , 

p = ±2u n{ r , _ ( 49 ) 

1 — re sin q = dn(t, re i ) at n>n c , 



where re 2 = Ti,/H c ; cn(t, re), dn(t, re x ) are the Jacobian elliptic functions shown in Fig. 3. 
With accounting for Eqs.(41), these expressions pass to Eqs.(37), (39) at H = H c . 

To elucidate the system behaviour with energy increase, let us expand the dependences (49) 
into Fourier series [2] 



a„cos[(2n — l)ut] at H < H c , 
p = ±8u{ n=1 x (50) 
I + J2 a-n cos(nut) at H > H c , 

n=l 

where one denotes 



at H>H C 



a n = i+g*- 1 d ; £ ^ 7' (51) 



fc = exp (^"^ > ^' = ^(Vl - re 2 ), K = K(R), 
parameter re is determined by Eq.(47). According to Eqs.(43), one has asymptotics 



, i re /32 at re 1, 

fc ~ U « exp(-7r/iV) at 1 - re 2 < 1, 
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where number N = uq/u is asimptotically as follows: 



N 



1 at H« H c , 

±ln^g^ at \H-H C \<^H C . 



(53) 



Thus, at low energies a single harmonics prevails to correspond to the coefficient of the 
Fourier series (50) 



On ~ 



H 



32H 



n-l/2 



H<:H C . 



(54) 



Respectively, near the critical energy, where the Fourier series gets the harmonics number 
N 3> 1, one obtaines 



a n ~ 8a; 



exp[— ir(n/N)] 



at 1 < n < N, 
at n> N. 



(55) 



As is known [2], above described behaviour is characterized by nonlinearity parameter 



a = 



dlnu; 



din/ 



(56) 



According to (42), (46), this parameter is determined by the following equation: 



a = 



1 E{k) 



l-K 2 K(k) 



1 



K 2 ( EjK- 1 ) 



i) \ 2 

), 



at k < 1, 
at k > 1. 



(57) 



The nonlinierity parameter takes on value a = at = and then arises indefinitely at the 
critical energy H c , tending to magnitude a = 1 at H — > oo. Such a behaviour is characterized 
by the following asimptotics: 



4 0-£r>T 



j(W/«c) 



16 



-H/W c 
16 



-2 r 



-2 



1 _ 1-H/Hc 1n 16 
1 2 111 1-H/Hc 



1 + 



H/Hc-1 



H/Hc-1 

1 - W c /W 



In 



16 



H/Hc-1 



at ?f« W c , 
at < 7i c - 7i < 7i c , 

at < 7i - 7i c < 7i c , 
at 7i > W c . 



(58) 



As a result, the observed picture of nonlinear oscillation is as follows. At low energy, when 
H <C H c , only single harmonic with the frequency uj ~ cjq keeps in the Fourier series (50), 
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so that the low-energy limit reduces to above considered case of harmonic oscillation (see 
Subsection 3.1). With energy increase, the harmonics number N arises in a manner of the 
dependence K(k) shown in Fig. 3, taking logarithmically large magnitudes (53) near the 
critical value Ti c = 2oJq. On the other hand, the oscillation frequancy (48) and the harmonic 
amplitudes (55) decrease monotonously to zero. Thus, one can mean in coarse manner that, 
with energy increase in the domain < 7i < Ti c , the single harmonic oscillation transforms 
to a harmonics superposition, whose number N increases monotonously to infinity, whereas 
frequency u> and amplitude a n ~ 8u decrease to zero. Just under the critical energy (0 < 
H c — H <C H c ) the system behaviour is characterizes by a set of solitons of different signs, 
whereas just above H c the signs of these solitons become equal (see Eqs.(49)). Remarcable 
peculiarity of such solitons set is that the width of a single soliton is reduced to At ~ u;^ 1 , 
whereas a distance between them is NAt ~ u^ 1 ^> ujq 1 . For just critical energy TC C , the 
system behaviour is caracterized by the separatrix solution (39) that reduces to the single 
soliton. 



3.3 Double well potential 



Taking into account use in the stochastic resonance problem [1], let us consider now the 
oscillation in double well potential. We will show that, in terms of the variables action-angle, 
the difference with above studied case of the nonlinear pendulum is quite quantitive, but not 
qualititive. 



The basic model is presented by the potential 



F = ^{l-x 2 )\ (59) 

counting off the energy at equilibrium positions x — ±1, uq is bare frequency. In this case, 
Hamilton equations (22) take the forms type of Eqs.(35): 



x — p, p = uJqx(1 — x 2 ) . (60) 

Near the saddle point x = 0, p = 0, corresponding phase portrait has the form differring from 
the same for nonlinear pendulum. This form is characterized by the separatrix (cf. Eq.(37)) 



p = ±u x y/1 - x 2 /2, (61) 

corresponding to the condition H = (u /2) 2 . Similarly to the case of nonlinear pendulum, 
the first of Eqs.(60) arrives at the separatrix law of motion: 
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x = ±2 1/2 y / l - [sinh(a;ot)] 2 . (62) 
For the time-dependence of the momentum, one has double-soliton solution (cf. Eq.(39)) 



p = T 2 1/2 cj sinh(cu t)^/l - [sinh(w t)] 2 , (63) 

where choice of signs corresponds to solitons moving to right or left sides. According to the 
motion laws (62), (63) - on the one hand, and Eqs.(38), (39) - on the other, the difference 
between the separatrix solutions for double well potential and nonlinear pendulum is that, 
in the course of the time near the point t — 0, the momentum gains two peaks of different 
signs in the former case and the single peak in the latter. 

To introduce the variables action-angle, it is convenient to use parameter re type of given by 
Eq.(40) and new variable £ determined by equality type of (41): 



2 , f re sin £ at re < 1, , nr x 

x = 1 - < . (65) 
[ sm( at re > 1 . 

Moreover, we shall need in using integrals 



l n (a, k)= [ A— = = -d£, a < 1, n = 0, ±1, ±2, . . . at re < 1; 

•> v 1 — ? 
o 



? (1 _ K -2£2\i-ra 

J n (a,K)= df, a< 1, n = 0,1,2,... at re > 1. 

* v 1 - ? 



(66) 



The form of dependencies Z„(re) = X„(o; = 1, re), J n (n) = J n (a = 1, re) is depicted in Fig. 6, 
the corresponding asymptotics read: 



7T 1 

J «(«0 ~ ^ - - !)«> «<1; 



J"„(re) « 2 + — (2n- l)re~ 2 , re > 1. 
The integrals are subjected to simple derivation rules 



(67) 
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K lt = (2n " 1} (Jn ~ Jn - 1} ' 



(68) 



where the arguments a, k are supressed for brevity. 

As a result, the last of definitions (24) arrives at the expression (cf. Eq.(42)) 



r-(l-/t 2 )J + 2J 1 - J 2 at «<1, 
2V2 I at k>1. [W) 

The action I = (2/n)S follows from this at a = 1. Respectively, for the proper frequency 
one obtains instead of Eq.(46) 



R 2 / [-(1-k 2 )J_i + (3k 2 -1)2o + 5Ji-3^] at « < 1, 

w = ^ W 1 (2^) 1 at « > 1. ^ 

Energy dependencies 1(H), ct'(W), following from Eqs.(69), (70), are depicted in Fig. 7. It is 
seen that these take the form type of the corresponding dependences for nonlinear pendu- 
lum (see Fig. 4). Fourier spectrum of the time dependence p(t) of the momentum behaves 
in analogous manner as in Eqs.(50): with energy increase in the domain < 7i < H c , sin- 
gle harmonic oscillation transforms to a harmonic superposition, whose number increases 
monotonously to infinity, whereas frequencies and amplitudes decrease to zero. Such a be- 
haviour is characterized by nonlinearity parameter (56), taking the following form (cf. Eq. 
(57)) 



[(1-k 2 )J-2-2J-i+k 2 J +2Ji-J 2 ] [-(1-k 2 )J +2Ti-T 2 ] 

a = { [-(i- K 2 )x_i+(3 K 2 -i)x +5Xi-3X 2 ]" K ~ ' (71) 

!l Zfjfr at k > 1. 

Respectively, double curvature A = 2^y of the dependence H(I) is connected with the 
parameter a as A = 2(u / I)a to read 



I Rk 2 ( 1 ~ k2 ) j -2-2J-i+« 2 Jo+2Ji-J 2 + < -I 

27r 2 J [-(l- K 2 )X_ 1 + (3 K 2 -l)J +5Xi-3J 2 ] a dl ' C72) 

^ at «>1. 
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3.4 Nonlinear pendulum under constant external field 



Before studying the effect of periodical external field, let us anounce main peculiatities of the 
constant perturbation following supersymmetry theory [6]. In this case, the potential energy 



F = 2u% sin(x/2) - 8x 



(73) 



is characterized, becides the proper frequency u>o, by a field strenth 8. Swithching such bias 
field arrives at the expression for flux j = (x)/t as follows: 



. _ 2nT nS 
3 ' ~ 1 T>{8) I T 7 ' 

where one introduces the integral 



(74) 



1(8) 



F(x) 1 8 
exp \ — f dx, sinx e = — , 



T 



(75) 



T is the temperature. In ergodic systems, the diffusion coefficient determined by equality 



((x-jt)*)=Dt 
takes the form 



(76) 



Z) = T— = 'v ch — 

c£ 1 T>{8) T' 



(77) 



However, non-ergodicity effects arrive at much more complicated form of the diffusion coef- 
ficient [6] 



, tt8 1 .5 J£ 
vrcth- + -arcsm-, + 



(78) 



where X is characteristic magnitude of the integral (75). Dependencies j(8), D(8) related to 
Eqs.(74), (78) are depicted in Fig. 8. 
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4 Nonlinear resonance conditions 



With accounting dissipation effects shifting the bare frequency u according to Eq.(28), the 
resonance condition reads 

m n w = n, w = \l[u{H)] 2 - 2- 2 , n=l,2,..., (79) 



where ft is the frequency of external signal, oj(7i) is the energy-dependent proper frequency 
determined by Eqs.(48), (46), m n is a resonance multiplicity being a rational number related 
to the natural one n. Formally, this condition means that we have to consider the phase 
portrait plane, which is revolved with angle velocity m n w. Then, the external addition of 
the Hamiltonian (23) is written in the form 

Hi = -Aj bo costf, d = (ft - m n w) t, (80) 

where the internal power f 00 corresponds to resonance conditions (79). Respectively, the zero 
term of the Hamiltonian (23) can be expanded near the resonant action I as follows: 

Ho-Hoo + ujoo 5I + ^(5I) 2 , 51 = I — I . (81) 
Here resonant action J , as well as other resonant values 



= 2-^a , a = a(I = I ) (82) 
Jo 

are given by resonance condition (79). It is worthwhile to note that curvature A/2 is deter- 
mined by nonlinearity parameter (56) taken at resonance condition I = Iq (see Figs. 5). The 
energy-dependence of the curvature is characterized by the following asymptotics: 



T^oo = 7~(-o(I — Iq)i ^oo = 



dfto 
dl 



An EE 2 



7=/n 



duj 
dl 



(27T)" 



1 — 12L 



n 2 < 



H/Hc-l 

2tt 2 V 1 H ) 



1 _ 1-H/Hc i n 16 
x 4 111 1-H/Hc 

1 ^ 4 111 n/Hc-i 



at H« W c , 

at T~L C — Ti <^ Ti c , 

at H - H c < W c , 
at 7i > W c . 



(83) 



As a result, Hamilton equations (25) arrive at the equation type of Eq.(36) for nonliniar 
pendulum: 

# + c£sin# = 0. (84) 
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Here, instead of the proper frequency uq, value u m stands for modulation frequency being 
determined as follows: 



oo 



\ 1/2 



V 2 / 



(85) 



Respectively, maximum value of the resonant energy variation and corresponding magnitude 
for the action 

5H m = Af 00 , 6I m = (^] ' (86) 



are determined by Eqs.(80), (81) to characterize a resonance window. 

Thus, we obtain the following picture of nonlinear resonance. At given magnitudes of the 
frequency f2 of external signal, the condition (79) fixes the system energy H as follows: 

71 tw \ f 1 at k < 1, 

— = K{k) { _! - ' 87 

Uo 2 +(2n/m n u o y l« at k>1. 



According to Fig. 9, the corresponding dependence 7ioo(^) has two branches, the lower of 
which relates to finite motion, the upper - to infinite one. Energy 7ioo related to the former 
falls down monotonously from the upper magnitude fixed by condition K(k) = hujq to zero 
within interval < Q < Q m , where maximum frequency is 

Vt m = m n u} \Jl - (2uj )- 2 . (88) 



Respectively, energy of the infinite motion arises monotonously from minimal magnitude 
fixed at Q = by condition k^ 1 K(k^ 1 ) = nu to 7i o ^ oo at ^ oo. Using abbreviated 
notation 

Q 2 = 1 + {2/m n ) 2 Q 2 , (89) 



one obtains the following asymptotics: 



H 



oo 



1 — 16exp — 2n(uj /Cl) 
4[l-(2^r 2 ]^ 



at < Q m — n <C Q m 



(90) 



for finite motion and 
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H 



oo 



He 



l + 16exp — 2n(iu /il) at Ocfi m , 



(_Q_Y 



at n > VL r 



(91) 



for infinite one. Frequency-dependencies of the resonant magnitudes Iq(Q), u>oo(Q), a(Q), 
Aq(Q) of the action, the proper frequency, the nonlinearity parameter and the double cur- 
vature of curve Hq(I) are depicted in Figs. 10, 11. In the case of the finite motion, these 
dependencies are characterized by the following asymptotics: 



f «5 OLUq 



n- 1 {l - 87r(cj /fi) exp -2-k(ojq/SI) } at Q < ft m , 



[1 - (2w ) 



— 2 ] S^ttT. ^ 



at < Q m — Vt -C fi m ; 



(92) 



^00 ~ ^0 



at Q <^ Q r 



n 

1 - [1 - (2u; ) -2 ] %^ at 0<fi m -fi<fi m ; 



(93) 



A 



i - [i - (2cj )- 2 ] at o < n m - n < n r , 



Respectively, for the infinite motion one has: 



(94) 



I pa 2uj 



| {l + 87r(cj /^) exp -27r(cj /^) } at Q < fi m , 



(95) 



cjqo ~ ^/2 at f2 <§; fi m and f2 3> f2 m ; 



(96) 



Ao 



1 

2 



l fjM 2 / 



2 \ 2u)q / I 87rwo 



exp 



(2a; /fi) 1 



27r(cj /fi) + l} at J) < Q m , 



at Q > fL 



(97) 



At u ^> 1/2, when Q m pa m n u , effective frequency fl may be replaced by the short-cut 
term (2/m n )Vt S> 1. 

If the oscillations were linear in character, the resonance might realize within window located 
between energies 7ioo — 57i m and 7ioo + In this window, the oscillations would have 

single frequency w determined by Eq.(79) and would be modulated with frequency oj m given 
by Eq.(85). However, nonlinearity effects described in Subsection 3.2 arrive at expanding 
the harmonics number to magnitude N = u /lj oq > 1 and narrowing the energy window 
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to width Hm = 2cj^. So, nonlinear resonance is realized for a share of nonlinear oscillations 
determined by the ratio 

— — = A = 2— a , (98) 



where Eqs.(82), (85), (86) are taken into account. The double curvature A(f2) decreases 
monotonously with the external frequency growth within the interval < Q < pQ m . This 
means that preference of the nonlinear resonance decreases with this frequency growth. 



5 Stochastic resonance conditions 



The stochastic resonance is known to be observed at condition (79), where m n = (2n)~ 1 , 

n = 1, 2, This condition means that during a period T = 2it/Vt of the external oscillation 

the stochastic system has a time to overcome an energy barrier AF even times 2n [1]. The 
proper frequency of stochastic resonance uj = 2tttk is reduced to the Kramers' rate rx given 
by Eq.(2). This case differs from the above considered case of nonlinear pendulum by only 
inserting temperature T instead of the system energy TC. In the case of double well potential 
considered in Subsection 3.3, the energy barrier AF = (tu /2) 2 is related to the temperature 
by characteristic parameter k 2 = 4T/u 2 . Then, the frequency of the barier overcoming is 
defined by equation 

lu = u K exp(-/sT 2 ), k 2 = AT/lu 2 (99) 



instead of the corresponding equation (46) for nonlinear pendulum. Taking into consideration 
dissipation effects, we obtaine the following condition of the stochastic resonance: 



T n = 



In 



1 + (4n) 2 n 2 



(100) 



As is shown in Fig. 12, the resonant temperature T (fl) arises monotonously, taking indefinite 
values at characteristic frequency 

n m = - (2u K )~ 2 . (ioi) 



Basing on the dependence T (Q) replacing above used relation 7i(f2), we are in position to 
consider the stochastic resonance conditions in analogy with nonlinear pendulum. 

Here, instead of Eqs.(40), (42), (46), the equalities (64), (69), (70) determine the resonant 
magnitudes Iq, uoq, a, A of the action, the proper frequency, the nonlinearity parameter 
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and the double curvature of the resonant dependence 7ioo(I)- Corresponding dependencies 
on the temperature are depicted in Figs. 13, 14 to show indefinite increase of the proper 
frequency and the action with tending to the characteristic magnitude (101). Related values 
of the nonlinearity parameter and the curvature decrease thereby. 



6 Conclusion 

The observed picture of nonlinear oscillations shows that with energy increase the sin- 
gle harmonic oscillation transforms to a harmonics superposition, whose number increases 
monotonously to infinity, whereas frequency and amplitude decrease to zero. In other words, 
with tending to a threshold energy, transition of the harmonic oscillations into a set of soli- 
tons is observed. We have shown that such behaviour, being typical for both the nonlinear 
pendulum and the double well potential, is characterized by nonlinearity parameter (56). 

Our choice of the action S and angle d as principle variables is caused by that Hamiltonian 
of free nonlinear oscillations does not depend on the angle d. In accordance with Eq.(84), 
switching on external harmonic signal arrives at a modulation of nonlinear oscillations with 
characteristic frequency (85) and nonlinearity parameter (82). It is appeared that nonlinear- 
ity effect narrows the resonance window to width fixed by the ratio (98) that is reduced to 
the curvature A of the dependence H.(I). According to the dependence Ao(fi) depicted in 
Fig. 14, this window is shrunken with the external frequency growth. As a result, preference 
of the stochastic resonance decreases with growth of the frequency fl of external signal. 
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Figure captions 

Fig. 2. (a) Motion laws for nonlinear pendulum (curve 1) and double well potential (curve 
2); (b) Corresponding time dependencies of the conjugate momentum. 

Fig. 3. Form of Jacobian elliptic functions and integrals. 

Fig. 4. Energy dependencies of the action and the proper frequency for nonlinear pendulum. 
Fig. 6. Form of the integrals (66). 

Fig. 7. Energy dependencies of the action (a) and the proper frequency (b) for double well 
potential. 

Fig. 8. Dependence of the flux j and the diffusion coefficient D on the force E of external 
field for nonlinear pendulum. 

Fig. 9. Dependence of the nonlinear pendulum energy on the external signal frequency: curve 
1 relates to finite motion, curve 2 - to infinite one. 

Fig. 10. Frequency dependencies of the action (a) and the proper frequency (b) for nonlinear 
pendulum: curve 1 relates to finite motion, curve 2 - to infinite one. 

Fig. 12. Dependence of the stochastic resonance temperature on the external signal frequency. 

Fig. 13. Frequency dependencies of the action (a) and the proper frequency (b) of the stochas- 
tic resonance. 

Fig. 14. Frequency dependencies of the nonlinearity parameter a and curvature A of the 
dependence 'H(I) at stochastic resonance. 
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